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ABSTRACT 


Let R be a regular local ring and take w to be an isolated singularity on R. Taking 
Z/2-graded R-modules, X and Y, a matrix factorization of w is a pair of morphisms 
(y, w) such that you = w and wow = w are satisfied in the diagram X > Y ee 
We will discuss the category of matrix factorizations of w in R and lead into the 
homotopy category of matrix factorizations as well as its historical development. 
Finally, we will conclude with the statement of the Kapustin-Li formula for the duality 
pairing on the morphisms in the matrix factorization category of (R,w) and discuss 


its implementation in SageMath. 
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CHAPTER 1 


INTRODUCTION 


Matrix factorizations were introduced by Eisenbud. He showed in [6] that by taking 
a finitely generated maximal Cohen-Macaulay module over the ring R/(f), where R 
is regular and local and (f) is a principal ideal in R, its minimal free resolution is 
obtained from a matrix factorization of f in R. 

The homotopy category of matrix factorizations of f in R was established by Aus- 
lander in [2] as a Calabi-Yau category which gave rise to the following interpretation 
by Kapustin and Li: considering the homotopy category, which is triangulated, as the 
category of boundary conditions in the Landau-Ginzburg B-model that corresponds 
to (R, f) allowed the derivation of a formula for the duality pairing on the morphism 
complexes in the matrix factorization category of (R, f), discussed in [11]. 

The focus of this thesis is the implementation of the Kapustin-Li Formula in an 
affine setting. In order to understand the statement and the code related to it, we 
start by introducing basic algebraic and complex analytic definitions and theorems 
in Chapter 2. We also look at introductory category theory and begin the intuition 
for the homotopy category. Then, in Chapter 3, we discuss matrix factorizations, 
the category comprised of them as objects, and the homotopy category of matrix 
factorizations. We end the chapter with the statement and a proof of Eisenbud’s 
matrix factorization theorem from [13]. 

We then move to the statement of the formula and its interpretation in Chapter 
4. Finally, discussion of the implementation of the Kapustin-Li formula in SageMath 


through the writing of Python scripts is our concluding chapter of content. The 


4 case can be seen in full in Appendix 


computation in the n = 1 variable of the form x 
A and is discussed in Chapter 5. For n 4 1, a complete calculation of the pairing could 
not be demonstrated, though discussed and attempted. The final chapter expresses 


the desire for future work concerning the formula to include the full implementation 


for any n-variable singularity. 


CHAPTER 2 


DEFINITIONS AND PRELIMINARY KNOWLEDGE 


First to be introduced will be key definitions from algebra, category theory, and 
complex analysis. Formulations of these definitions can also be found in the standard 
references ([3], [5], [10], [1]). As they arise, important and relevant theorems and 


propositions will be presented. 


Notation 2.1. The identity element of a ring R or R-module M will be denoted idz 


or idyy, respectively. 


Notation 2.2. The Krull dimension of a ring R we denote as dim R. 


2.1 ALGEBRA 


Remark 2.1.1. A common statement on the Krull dimension is the following of 
dim R = 0 if and only if every prime ideal P in R is a maximal ideal. A notion which 


follows is: if R is Noetherian, then FR has finite length if and only if dim R= 0. 


Definition 2.1.2. Let R be a Noetherian local ring with maximal ideal 
m = (Q1,...,@)) where n is minimal. Then R is regular if dim R =n. Here we refer 


tO a4, ...,dn as a regular system of parameters. 


Definition 2.1.3. Given additive subgroups R; C R such that R = @, R; and 
RR; C Riz;, then we say R is Z-graded, or Z/2Z-graded. 


This property of gradedness can also extend to the modules of the ring under 


certain conditions. 


Definition 2.1.4. Let R be a graded ring. Given additive subgroups M; C M such 
that M = @; M; and R:M; C M;i+;, then we say M is Z-graded over R. 


Definition 2.1.5. For a R-module M, a sequence of elements ay,...,a, is a regular 
sequence on M if (a1,...,@n)M #4 M and, for i = 1,...,n, a; is a nonzero divisor on 


M/(ay, .--; @n—1)- 


We have seen the geometric measure of a module in the idea of Krull dimension; 


now we may look at the homological measure of the size of a module M. 


Definition 2.1.6. If R is local Noetherian with maximal ideal m and M is a nonzero 
finitely generated R-module, then grade(m, /), or the depth of M, is the length of 


the longest regular sequence in m on M. 


Definition 2.1.7. A R-module M is a (maximal) Cohen-Macaulay module if the 


depth of M is equal to the Krull dimension of R. 


We denote a Cohen-Macaulay module by the abbreviation CM. The ring R is a 
Cohen-Macaulay ring if Risa CM over R. 


Definition 2.1.8. An exact sequence is a sequence of objects F; and morphisms 
between these objects f;, written as 

| fae a CE eR 
such that the image of a morphism is the kernel of the following morphism, 


im f; = ker fi41. 


Definition 2.1.9. A resolution is an exact sequence of modules. A free resolution 


is one where each module F; is free. 


Definition 2.1.10. Given an exact sequence of R-modules, 
03 N > Fh-1 O.... 9 F, OM ROM 0 


for F; free R-modules, then N is the n-th syzygy of M. 


Definition 2.1.11. An injective resolution of the module M is an exact sequence of 
the form 


0+ M3 Pore... (21) 


where J’ are injective modules. 


An important conclusion is that every module M has an injective resolution. The 
dual notion of these resolutions is the projective resolution. Every module N also has 


a projective resolution, an exact sequence of the form 
. 273 Ph HRN (2.2) 


where P; are projective modules. Definitions of injective and projective modules can 


be found in [3] also containing propositions, with proof, relaying their properties. 


Definition 2.1.12. For a module M which admits a finite injective resolution, the 
minimal length among all finite injective resolutions of M as seen in (2.1) is its 


injective dimension, denoted injp(M). 


Definition 2.1.13. Similarly, if a module N admits a finite projective resolution, 
the minimal length of all finite projective resolutions of N as seen in (2.2) is its 


projective dimension, denoted pdp(V). 


Considering these homological measures, we come to a formulation by Auslander- 
Buchsbaum. The equation below can be seen in [5], [6], and [13] and shown in [12] 
where first two preliminary lemmata are proven and then used in the proof of the 


statement. 


Theorem 2.1.14. (Auslander-Buchsbaum) For R a commutative Noetherian local 


ring and M a nonzero finitely generated R-module of finite projective dimension, then 


pdpM = depth R— depth M (2.3) 


Definition 2.1.15. A Gorenstein ring is a commutative Noetherian ring such that 
each localization at a prime ideal is a Gorenstein local ring — a commutative Noethe- 


rian local ring R with finite injective dimension as a R-module. 


The following, which can be found in [13], are propositions detailing some useful 


properties of CM modules over a ring R. 


Proposition 2.1.16. Let R be a commutative Noetherian local ring and let 
0 — L — M — N — 0 be an exact sequence of R-modules. Then the following 


are true: 
1. If Land N are CM, then M isCM. 


2. If M and N are C'M, then so is L. 


Proposition 2.1.17. If R is a regular local ring, then any CM module over R is a 


free module. 


Finally, an important algebraic structure required for all proceeding material is 
that of the chain complex. There is a dual notion, cochain complex, however we will 


use the convention of chain complex for this paper. 


Definition 2.1.18. A chain complex (A.,d.) is a sequence of modules ..., Ag, Aj, ... 
connected by homomorphisms, called differentials, dy: An — An—1 such that 


dy, © dn, = 0. This can be represented by 
ge he Say ee As a (2.4) 
An important complex used in homological algebra is the Koszul complez. 


Definition 2.1.19. For R a commutative ring, EF a free R-module with finite rank 


r, and R-linear map y: E > R, the Koszul complex associated to vy is the chain 


complex of R-modules 
Oe he A oe 0 (2.5) 


where d; is the differential between the exterior powers of the free module; note, 
A°E = R. For any e; € E, dz is defined as 
k 


di(er A... A ex) = S(-1)*tHe A... A ples) A. A ee (2.6) 
i=l 


2.2. CATEGORY THEORY 
In this section, we introduce basic category theory, including additive categories and 


triangulated categories. These concepts are important to and applied throughout the 


contents of Chapters 3 and 4, and each reference therein. 


Definition 2.2.1. A category @ consists of: 
1. A class of objects Obj(@), 


2. A class Homy of morphisms between two objects, where each morphism has a 


source and a target in Obj(@), and 


3. For three objects A, B,and C, a binary operation called the composition of 


morphisms where Hom¢(A, B) x Hom¢g(B,C) > Hom¢(A,C), f x go gof. 


such that there is exactly one identity morphism for every object and that associa- 


tivity of morphisms holds. 


Definition 2.2.2. Let A and B be two objects in a category @. A product of A and 
B is an object P along with morphisms A “~ P —*» B such that, given any diagram 


A X —) B, there exists a unique morphism u : X —> P so that the following 


diagram commutes. 


xX 
ae Ne 
oa 7 P2 ee 


Definition 2.2.3. Let A and B be two objects in a category @. A coproduct of A 


(2.7) 


and B is an object C along with morphisms A 4 C <  B such that, given any 


diagram A “> Y  B, there exists a unique morphism u : C —> Y so that the 


B 
©2 (2.8) 


following diagram commutes. 


Y1 


A 


Y 
u 
Cy C 
A coproduct is the dual notion of a product, and they are often denoted AJ|] B 
and AJ[[ 8B for products and coproducts, respectively. 
Just as in ring theory, kernels and cokernels also have applications in category 
theory, though the added context of an object is necessary to understand more than 


just how the morphism acts. There is also an assumption that the category contains 


zero morphisms. 


Definition 2.2.4. For a category @, let a: X — Y be some morphism between 
objects X,Y in @. Then for any morphisms g,h: A — X for some object A in @, if 


ag = ah we call a the zero morphism. 


Definition 2.2.5. For a category @ which contains zero morphisms and for some 
morphism f : X — Y, the kernel of f is an object P in @ with the morphism 
p: P — X, written ker f = (P,p), such that the composition f o p is the zero 


morphism from P to Y, which we denote Op in the following diagram. Explicitly, 


given any p’: P’ + X such that f op’ is the zero morphism, there is a unique 


morphism u: P’ + P such that pou =p’. 


(2.9) 


Definition 2.2.6. Similarly, for @ which again contains zero morphisms and some 
morphism f : X — Y, the cokernel of f is an object Q in @ with the morphism 
q:Y — Q, written coker f = (Q,q), such that the composition qo f is the zero 
morphism from X to Q, denoted Og in the following diagram. Explicitly, given any 
qd: Y + Q such that q' o f is the zero morphism, there is a unique morphism 


u:Q—>Q such that uog=q'. 


(2.10) 


Definition 2.2.7. A category @ is an additive category if the following hold: 


1. For every X,Y € Obj(@), Hom(X, Y) is an abelian group and the composition 


of morphisms is bilinear, 
2. @ contains a zero object (an object that is both initial and terminal), and 


3. For any X,Y € Obj(@), there exists a coproduct X [[ Y in @. 


A functor F : @ + JY between two such categories is additive if for all X,Y € Obj(@) 


F induces a homomorphism of groups Homy(X,Y) + Homg(X,Y). 


Definition 2.2.8. Let 7 bean additive category and let % : FY > Y bean additive 
automorphism. A triangle in Z is a sequence A, “> Ay “> A; “*> XA, of objects 


and morphisms in 7. 


Leb As Ae ES A and By ean Bo see, B3 Beck “B, be two triangles in 


ZT. A morphism of triangles is a commutative diagram. 


a a a 
Aj : > Ao > A3 z me MA, 
Y1 p2 3 UY] 
B, By 2 By Be . Bs Bs . SB, 


(2.11) 


If ~1, 2, andy3 are isomorphisms in 7, then (1, 2, Y3) is called an isomorphism of 


triangles. 


Definition 2.2.9. Let Y be an additive category. Then 7, an additive automor- 
phism %, and a collection A of distinguished triangles is a triangulated category if all 


of the following are satisfied: 
1. If a triangle is isomorphic to a triangle in A, then it is in A, 
2. For every A € Obj(.7Z) the triangle A —> A —> 0 — DAisin A, 


3. For every A,, Ag € Obj(7) and a € Homz(Aj, A2) there is a triangle in A of 


the form Ay a A» =>? A3 ==>? wAj, 


Ae Bor Aj Ay 8 Ay 8 SA ii As then Ay 5 Ay =) SA DAS is 


in A, 
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5. Given two triangles A; “> A, > As 8» SA, and 
B, as By ae Bs aes SB, in A, each commutative diagram can be completed 


to morphisms of triangles. 


a a a 
Aj : > Ao 2 > A3 : > MA, 
Y1 p2 3 UY] 
B, py ' By Be : Bs Bs ¢ EB, 


(2.12) 


Definition 2.2.10. Let @ be an additive category. Let C, be a chain complex with 
boundary maps de» : Cy —> Cp-1. For any k € Z, the k-shifted chain complet Ck]. 


is defined by doyx,n = (-1)*don+k as Clk]n = Cnn, written 

dejan: Clkln —> Clk] n—1 (2.13) 
Definition 2.2.11. Given two chain complexes A and B, and two chain maps, 
f,g: A— B,a chain homotopy is a sequence of homomorphisms hy: An > Brit 


so that 
f-—g=hda+dgh (2.14) 


This can be represented by the following diagram. 


dan dawn da,n dag 
we eee n+1 atlas = A, - > An-1 ie : eae 
Ans Fina | a | ie Tie | bg 
Bay ig B, > Bh-1 fos. 
dp n+2 dpn+i dgin dpn—1 (2.15) 


The map hd,4+d,h induces the zero map on homology for any h, thus f and g induce 


the same map on homology. The maps f and g are said to be (chain) homotopic. 


Definition 2.2.12. The homotopic maps define an equivalence relation on the 
abelian groups of morphisms in the category @ which we call the equivalence class of 


the morphisms. 
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Definition 2.2.13. Given an additive category of chain complexes, @, the homotopy 
category, denoted H(@), retains the same objects as the category @, however the 


morphisms are the equivalence classes of chain maps. 


2.3. COMPLEX ANALYSIS 


In this final section of definitions, we state concepts in complex analysis that will be 


used in Chapters 4 and 5. 


Definition 2.3.1. Let f be holomorphic everywhere except at a point zo. We say 


zo is an isolated singularity. 
Definition 2.3.2. The residue of f at zo is defined as 
1 
Res, f = a f f2)dz (2.16) 
271 


In calculating the residue, we can consider the series expansion of f about Zp, 


written: 


12> 3s Cn(Z — 29)”, where Cp = : ¢ ACS (2.17) 


n=—0o 2n0 (¢ aa a 


Thus Res,,f is the coefficient of (z — z9)~! in the expansion, so 


Rey feese = facs(o (2.18) 


This alternate residue calculation will make the implementation of the Kapustin- 


Li formula described in Chapter 4 easier and more explicit to calculate in SageMath. 
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CHAPTER 3 


MATRIX FACTORIZATIONS 


In this chapter, we will define matrix factorizations as well as the category which 
contains them as objects. Next, understanding of the homotopy category of matrix 
factorizations will be developed in order to understand our final portion of the chapter 
on the first contribution of matrix factorizations which can be referenced in [4], [7], 


[9], [13]. 
Definition 3.1.1. Let R be a commutative ring, with w € R. A matrix factorization 
(A, B,y, wv), or shortened (y, wv), of w in R is a diagram 


A B (3.1) 


for A, B finitely generated free R-modules and vy, ~ R-homomorphisms such that the 


following are satisfied: 
pow=w-idg and pop=w-idy (3:2) 


Example 3.1.2. Let R = C[z] and w = x”. Considering R as the R-modules, we 


have the factorizations 


where is just multiplication by x? and w by x”~¢ so that qo y = 2”. 
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Remark 3.1.3. We can see for matrix factorization (y,w) of w that w annihilates 
coker y as defined in (2.10). So w(coker y) = 0. This will be used in (3.1.4) to clarify 


a condition. 
The following propositions are shown in [6] 


Proposition 3.1.4. Let R be a Noetherian ring and ... > A Bs ees en 


be a free resolution of finitely generated R-modules which is periodic of period 2. 


Then rank A = rank B. 


Proposition 3.1.5. Let w € R be a nonzero divisor, and let y : A — B be a map 
between free modules. There exists a matrix factorization of the form (y, w) if and 


only if 
1. rank A = rank B 
2. det y = rank B, and 


3. w- Fit 1(y) C (det vy) 


Remark 3.1.6. The third condition above refers to the fitting invariant of coker y, 
denoted Fit;(y). However, Eisenbud discusses that the assumption of the proposition 
along with the first condition of equal rank implies the annihilator of coker y, written 
annp(coker y), is equivalent to (Fit;(y) : det y). So for w- Fit 1(y) C (det y) we 
may instead write w(coker y) = 0. The proof found in [6] uses this idea of condition 


three. 


Definition 3.1.7. A morphism 6 between two matrix factorizations (A, By, Y1, v1) 


and (Ag, Bo, Y2, 2) of w is a pair of maps a: Ay — Ag and 6: By — By such that 
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the following diagram commutes. 


fe a 
a B a 
Ap 2 - By CD) ae 


(3.4) 


The commutativity of the left side implies the commutativity of the right side 
of (3.4) so we can also state, though redundant, that a ow, = W206 must also be 


satisfied by a and 6 to be a morphism between matrix factorizations. 


Definition 3.1.8. Two matrix factorizations of w are equivalent if a and £ are iso- 


morphisms. 


Definition 3.1.9. Equivalent matrix factorizations which have non-unit maps are 


referred to as reduced matrix factorizations. 


For R a regular local ring, any matrix factorization can be written, using the 
differential (y,~), as the following direct sum with a reduced matrix factorization 


(yr, wy) and a,b € Z2° 
(y, b) = (Yr, Ur) ® (id 4, w)* ® (w, idg)? a5 


Specifically, if the matrix y contains a unit, then we can write it as the sum of a 
reduced matrix factorization and (id4,w) to some power a. Similarly, if ~ contains 
a unit, we can write it as the sum of a reduced matrix factorization and (w,idg) to 


some power b. 


Example 3.1.10. Again considering the example in one variable of the matrix fac- 


d Re) 


torizations of x”, written (7°, x2"~°), we see for d 4 0 the matrix factorization of w is 


reduced. 
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Now we can build the category of matrix factorizations to have the collection 
of matrix factorizations as objects and to let the morphisms between the matrix 


factorizations as defined in (3.1.7) be the morphisms of the category. 


Definition 3.1.11. The category of matrix factorizations M F(R, w) is the collection 


of matrix factorizations of w in R and the morphisms between them. 


MF(R,w) can be observed to be an additive category with the expected zero 
object and direct sums as the coproduct. 

Next we introduce the notion that the morphisms between matrix factorizations 
can be chain homotopic and the conditions required to be such. This will be used in 
each remaining chapter and is important to the idea of the Kapustin-Li formula as 
well as the major theorem in this chapter. Each component of the morphisms between 


matrix factorizations will have to satisfy applications of (2.14) to be homotopic. 


Definition 3.1.12. Let 0,6’ : (Ai, Bi, 1,01) > (Ae, Bo, v2, 2) be two morphisms 
in MF(R,w) where 0 = (a, 3) and & = (a’, 8’). Then 6,6’ are homotopic if there 


exist maps s,t which satisfy 


a-a’=soyityyot (3.6) 
B-B =tom+pmos (3.7) 
seen in the following diagram. 
At 1 Le Vy me 
a’,a a, ra al, a 
Ap 2 op. We 7 e 


(3.8) 


These homotopic maps define an equivalence relation on the abelian groups of mor- 
phisms in the category M F'(R,w) and we denote the equivalence class of a morphism 


0 by [6]. 


16 


Next we introduce the homotopy category of matrix factorizations which has the 
collection of matrix factorizations as objects and those homotopic equivalence classes 


described above as morphisms. 


Definition 3.1.13. The homotopy category HMF(R,w) is the category which re- 
tains the same objects as M F(R, x), but the morphisms are the homotopy equivalence 


classes of morphisms. 


It can be observed that HMF(R,w) is an additive category. The morphisms, 
which are the homotopy equivalence classes, form an abelian group and the compo- 
sition is bilinear; these maps act as matrix multiplication which is bilinear. Finally, 
the zero object and coproduct are retained from the matrix factorization category. 

Next, we provide the concluding material necessary to understand Eisenbud’s 
theorem on matrix factorizations. 

We first define a quotient which will lead to the categories used in Eisenbud’s 


matrix factorization theorem (3.1.15). 


Definition 3.1.14. Let @ be a category with the homomorphism sets being abelian 
groups, and & be a set of objects in @. We can now define the category @/& as 
the category which retains the same objects but whose morphisms between objects 
A and B are the elements of the quotient Homy(A, B)/@#(A, B), where (A, B) are 


all morphisms from A to B which factor through direct sums of the objects of /. 


Letting ¢ represent a direct sum of objects in ./, we have the following diagram 


describing the behavior of the morphisms (A, B). 


A > B 


Le 


The zero objects in & are retained in the category @/@ which can be shown to 
be an additive category. 

For R a regular local ring and (f) a principal ideal in R, we let S = R/(f). Then 
we can define @ (5) to be the category of all CM modules over S and the quotient as 
€(S) = G(S)/{S™} for all n. @(9) has morphisms in the quotient Homg over a S- 
submodule of Homg. Explicitly, let @(M, N) be the set of S-homomorphisms of M to 
N which factor through a free module, F’, written as M > F > N. This A(M, N) is 
a submodule of S so the quotient is written Hom,(M, N) = Homs(M, N)/A(M, N). 

We also note that the quotient MF(R, f)/{(idpm, f)} for all n is denoted 
MF(R, f). Considering RM F(R, f), the reduced matrix factorization category, the 
quotient RM F(R, f)/{(idpm, f), (f,idgm)} is denoted RM F(R, f), where (idg, f) 
and (f,idp) are obviously not reduced matrix factorizations as discussed and seen in 
(3.5): 

Finally, we define the additive functor coker : MF(R, f) > @(S), where coker 
(a, 3) is a homomorphism of CM modules coker (v1, 1) — coker (y2, 2) and the 
object associated to coker (idr, f) is the zero object and the object associated to 
coker (f,idp) is S. 

Now we are ready to state the following theorem. There is a proof outlined in [13] 


which follows. 


Theorem 3.1.15. (Eisenbud’s Matrix Factorization Theorem) Suppose R is a reg- 
ular local ring and (f) is a principal ideal. If S = R/(f) is a hypersurface, then coker 
induces an equivalence: 


MF(R, f) = @(S) (3.10) 


Furthermore, RMF (R, f) = @(S). 


Proof. As coker (idr, f) = 0, coker induces the functor MF(R, f) — @(S), which 


will be denoted Coker. For a nontrivial CMM module M we have a free resolution of 
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0— RM +; R™ — M0 
and we can obtain (y,w) which satisfies conditions in (3.2) where pow = f - idpm) 
and Woy = idgm:f. In order to find (y, W), we first note that since M is a S-module 
it follows that fM = 0, so for any x € R™ there is a unique element y € R\) where 
f-x=(y). Letting y = v(y), then we see w is a linear mapping from R™ to itself 
and therefore satisfies p-w = f-idpim. This defines a functor F : @(S) > MEF(R, f). 
So we set F'(M) = (y,W), which is determined uniquely as an object in MF(R, f) 
since we may neglect (idr, f), and we note that if we make the choice for y to be 
minimal and if (y1, 71) is another matrix factorization obtained from M, then there 


are invertible matrices a and ( such that the following is a commutative diagram 


; 
0 Raut Rt > M > 0 
B a | 
y 
0 » pm —*! pm eg | sr, 


(3.11) 


with y = (% ae) and 6 = (es foe) Therefore (a, (@) is a morphism from (7¥, 6) 
to (y1,W1). Now given a morphism g : M, > Mp in @(S), there is a commutative 


diagram 


0 > Rm Se M, - 0 
B a | 
” 
0 i a, - 0 


(3.12) 


Therefore (a, 3) gives a morphism of functors F'(M,) > F(M2), denoted F'(g). 
If we take (a’, 6’) to be another morphism which makes (3.13) commute, then we 
can define a homotopy uz: R™! — R™ such that (3.7) and (3.8) are satisfied. So the 


morphism (a, 3) > (a’, 6’) is a composition (1, W- Y~1) o (Y2, idr), where 
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(1, P1) : (91, V1) > (idr, f-idg) and (Yo, idg) : (idr, f-idp) > (Ye, 2). Therefore 
the morphism is in MF(R,f). Thus F(g) is uniquely determined. Then we can 
check that F'- Coker = idg and Coker - F' = idg which shows the desired equiva- 


lence MF(R, f) = @(S). Similarly, since coker(f,idg) = S, the second equivalence 


RMF(R, f) = @(S) follows. 
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CHAPTER 4 


KAPUSTIN-LI FORMULA 


Here we will first introduce the statement of the formula discovered by Kapustin and 
Li for the duality pairing as stated and discussed in [11]. We then follow with an 


accessible description of the formula. 


4.1 STATEMENT OF FORMULA 


According to Kapustin and Li, the formula provides the duality pairing on the mor- 
phism complexes in the matrix factorization category of an isolated hypersurface 
singularity. 

We consider a regular local ring R with isolated singularity w found by taking a 
maximal ideal m in R with w € m. A matrix factorization (A, B,y,w) of w in R in 
this context consists of Z/2-graded finite free R-modules A and B equipped with an 
odd endomorphism d which satisfies d? = w. This object then corresponds to the pair 


of square matrices y and wW which can be combined into a supermatrix 


0 
Q = 
w 0 
where Q? = w- id. So we see the odd endomorphism d = y @ w = (y, 1). 
Furthermore, we take the homotopy category of matrix factorizations, denoted 
HMF(R,w). For X,Y matrix factorizations, HMF(R,w)(X,Y) denotes the mor- 


phisms, more explicitly the homotopy equivalence classes, in the homotopy category 


between X and Y which we have called a and ( for Q in (3.8). 
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Similarly, HMF(R,w)(Y, X[n]) denotes the morphisms in the homotopy cate- 
gory between Y and X|n]. X|n] denotes the shifted complex of the matrix factor- 
ization X as described in (2.13) and to be shown in the following diagram. For 
X = (X°,X!,ox,yx) and Y = (Y°,Y1,yy,yvy), and for the homotopy equiva- 
lence class [ex]; say (ax, bx), from HM F(R,w)(X,Y) and for [oy], (ay, By), from 
HMF(R,w)(Y, X|n]) we have 


¥ Ox sei Wx Ja 
ax Bx Ox 
yo ea by Be 
ay By ay 
n a nt] = yn 


(4.1) 


Then we have the following formula. 


Definition 4.1.1. For some F in HMF(R,w)(X,Y) and for some G in 
HMF(R,w)(Y, X[nJ), 

1 tr(FG(dQ)) 
(277)"n! — O1WOQW...OnW 


(FG) +> 


The following theorem discusses the property of the non-degeneracy. The discus- 


sion and proof can be seen in [8]. 


Theorem 4.1.2. The formula defined above is a non-degenerate pairing, satisfying 


1. For F € HMF(R,w)(X,Y), if 6(F,G) = 0 for all G € HMF(R,w)(Y, X[n]) 
then F = 0, and 


2. For G € HMF(R,w)(Y, X[n]), if (F,G) = 0 for all F € HMF(R,w)(X,Y) 
then G = 0. 


Ze 


Finally, letting k = C allows the application of path integral methods to find the 


pairing. 


4.2. UNDERSTANDING THE FORMULA 


In this section we will discuss each component of the formula to be encoded and 
detailed in the following chapter. 

Given the morphisms, F = (ax,(Gx) and G = (ay, Sy), as shown in (4.1), we 
form matrices with the maps as entries. The following are for the singularity in odd 
n variables, where G is shifted by n, and thus has a matrix representation with entries 


within the off diagonal. 


Then we form the matrix 'G 


FG= 
Bx © By 0 


Then we take the differential (yy, Wy) and form the matrix 


Ox 0 yy 

dy 0 
We use the differential associated to Y as it is common to both morphisms, F’ and 
G. Then the matrix Q is differentiated by taking partial derivatives of the entries in 
each of the n variables of w, denoted 0,w for the i” variable, and wedging the forms 


to construct (dQ)’", shown explicitly as 


(dQ) = Oy Ot (4.2) 


Then matrix representations of FG and (dQ)’” as detailed above are multiplied. 
With this new matrix, we take its trace. This sum will produce a multivariate poly- 
nomial which then will be divided by the product of partial derivatives of the original 
singularity, seen in the formula as O,w...0,w. 

Finally, we are left with the need to integrate. As shown in (2.17) and (2.18), we 
can take the series expansion of our quotient and find the coefficient of the term with 
degree -1. This will be the evaluation of the residue, therefore as described in (2.16) 


we can rewrite our formula as follows 


{|Ojw|=e} O,wOow...d,W n) OW, Oow, ..., Onw 


i f tr(FG(dQ)*”) 1 ("f) Res tr(FG(dQ)*”) 


(Q7i)"n! = 


| (4.3) 
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CHAPTER 5 


SAGEMATH CALCULATIONS 


In this chapter, we will discuss the encoding of the Kapustin-Li formula using Sage- 
Math. Included will be excerpts of SageMath input and output as well as contents 
of Python scripts utilized as functions or methods for parts of the formula, with full 
references found in the Appendix. 

The conditions of the formula require R be a regular local ring. For the entirety 
this chapter we let R be the ring of formal power series in n variables with complex 
coefficients. We can declare this ring in SageMath via the following the PowerSeries- 


Ring() command, with a two-variable example shown below. 
sage: R = PowerSeriesRing(CC, [x,y]) 


This command outputs a description of the ring declared including the adjoined 


variables, the type of ring, i.e. Power Series, Polynomial, etc., and the field. 


Multivariate Power Series Ring in x, y over Complex Field with 53 bits 


of precision 


In general, the n-variable ring can be crafted using a while loop starting at index 
1 as seen below. The loop creates notation of n-variables as indexed «x’s for simpler 
code and stores them in a string to be used in the declaration of the ring R. This is 
useful in that there is a simple input of the number of variables desired and an output 
of the ring with the corresponding number of variables. This string of variables is 


also convenient to have as it can be called on for later calculation. 
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while index <= n: 
str = str+"x_{},".format (index) 
index = index + 1 

str = str i-1] 


R = PowerSeriesRing(CC, str) 


Similar to the previous example, the output describes the ring with the variables 


and the field it is over. In the case of n = 4, we have the following output. 


Multivariate Power Series Ring in x_1, x_2, x_3, x_4 over Complex 


Field with 53 bits of precision 


Since the formula begins with the input of morphisms F' and G, we need the 
singularity w and the matrix factorization, Y, in order to construct the matrix Q. 
We can establish F and G by constructing matrices from their provided maps as seen 
in (4.1). For the n = 1 variable case, the computation of (dQ)/” = (dQ) = dQ is 
easy and requires only differentiation of the entries in the matrix with respect to the 
single variable, here considered x. At the end of this chapter we have an example 
of a full implementation of the formula for n = 1 of the form 2% and specifically for 
the singularity «*. For n 4 1, we will skip the computation and construction of the 
differential matrix (dQ)/” and return to discuss it in Chapter 6. Now we switch to 
discussion of the residue calculation, given that a polynomial is the output of the 
trace calculation. This will utilize the taylor() command, included in SageMath, to 
form the Taylor expansion and the residue computation. First, we must establish the 
denominator of our polynomial in order to expand the expression in a Taylor series. 

As seen in (4.3), the trace evaluation is divided by the product of the partial 
derivatives of our singularity. With these declared n variables we can compute the 
partial derivatives. The general case Python script contains the while loop which 


takes the string of variables, parsed by variable and referred to below as "vari", and 
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takes the partial derivative then the product of it and all preceding partials. The 
index again begins at 1 and the product of the partial derivatives denoted "der" is 


initially declared 1. 


vari = str.split(",") 

while index <= n: 
partial = sing.derivative(var (vari Lindex-1])) 
der = der*partial 


index = index+i 


More simply, the one variable case script contains the command to differentiate 
the single variable singularity. Now we are ready to begin finding the Taylor series 
expansion. Given the polynomial output of the trace calculation, and the computed 
product of the partial derivatives of the singularity, taking the quotient we can find 
the Taylor series expansion. For the polynomial output, called "poly", and for the 
product of partial derivatives found above, still called "der", we can compute the 
quotient, referred "quot", by simple division: quot = poly/der. 

Then, we can use the taylor () command which takes 4 arguments: the function, 
the variables, the singularity, and the desired degree to which the terms be printed. 
The function to be used is our quotient, "quot". We have already established a list of 
variables so we implement that list "vari" for each a variable of the list in the second 
spot of the command. For the third argument, we want to find the solution set to 
each partial derivative, as seen in the formula statement in (4.1.1). Each solution set 
will yield values for that variable where the residue will be computed and summed 
for each. We can call the set an array and increment by 1 to evaluate the residue 
considering each entry. For the final argument, by the definition given in (2.18) we 
would like to find the coefficient of the degree -1 term of the series expansion so we 


place -1 in the final spot of our command. So to find our 7th residue the command 


yas 


would be as follows. 
taylor (quot,varila],soln[i] ,-1) 


To retrieve the coefficient from this expression, we can convert the output to a 
string and parse the initial term. Converting each of these parsed terms back to 
variables allows us to take their sum and complete our calculation of the residues and 
thus the complex integral. Though not completed, these commands and points were 
the focus of the attempt to implement the formula from the trace calculation on. 

We close the computation chapter with an example for the singularity w = 2’. 
The Python script used can be found in Appendix A. We let X be the factorization 
with maps (3,2) and Y be the factorization with maps (2?,x”). We can write the 


differentials in terms of the maps associated to X and Y, respectively. 


Per the formula, the matrix Y will hold the role of Q in the computation. As the 
variable count is 1, the matrix factorization associated to the shifted complex of X [1] 


is ee 


Finally, we define F' and G, consisting of maps of the respective homotopy equivalence 
classes, denoted (az, 3,) and (ay, Gy) respectively as discussed in section 4.2. Here, 
we let F' be defined by (x?, 7) and G be (x, 1) so the composition of F and G, denoted 
FG in the formula, forms a matrix of the multiplication of the maps as entries on the 


off diagonal, as Q is defined. 
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We then multiply the matrix FG by the differentiated matrix dQ, which is differen- 


tiated with respect to x by entry, associated to the factorization maps of Y 
0 
FG xdQ= x — 


Next, the trace of the matrix found above is computed and we form the quotient 
comprised of the trace output and the differentiated singularity, then we simplify the 
expression. 
Qe*+277 x 1 
4c Ci: «Oe 
The Taylor expansion of this polynomial about x = 0 is the polynomial itself so we 


simply take the coefficient of the 4 term as our residue. Therefore, for the singularity 
x* and for homotopy equivalence classes F = (x?,2) and G = (2,1), we have that 


the computed residue 5. 
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CHAPTER 6 


CONCLUSION 


We have discussed matrix factorizations and their history, the category of matrix 
factorizations, and their homotopy category. We also focused on Eisenbud’s Matrix 
Factorization Theorem as well as saw and discussed a proof. This intuition lead to 
a discussion of the Kapustin-Li formula as well as an effort to implement utilizing 
SageMath. A valuable next move would be to complete the affine application of the 
formula and develop a library of Python scripts for public use. 

Neither the construction of the matrix representations called F’, G, and the dif- 
ferential matrix (dQ) nor the final component of residue calculation in the general 
n # 1 case were completed in this project, though discussed and attempted. In 
particular concerning the differential matrix, differential forms and the wedging of 
elements in SageMath requires background of their syntax concerning the context of 
manifolds which was not discussed in detail here. In the furthering of this project, this 
understanding would be crucial and greatly aid in the development and completion 
of the formula’s calculation in this context and perhaps others. 

After completion, it would also be interesting and helpful to compile a library 
of sample singularities and matrix factorizations for users to develop a closer un- 
derstanding of the inner workings of the program. An ever-growing library could 
help to consider other statements and contexts leading to further study and potential 
research in this area. The completion of the implementation in SageMath of each 
portion of the formula in the context of the formal power series ring over the complex 


numbers would allow for quicker computation of the duality pairing and supply a 
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referable material for other contexts. 

Further research and implementation into the homological and topological under- 
standing of the Kapustin-Li formula, and by extension matrix factorizations and their 
homotopy category, could be taken in many directions. Hopefully as questions are 
raised and solved, the impact and understanding of the formula and its applications 


across mathematics and other sciences will expand. 
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APPENDIX A 


ONE VARIABLE 2” CASE 


The contents of the Python script which computes the duality pairing for any one 


variable case is below. 


from sage.all import * 


## Here we assume a single variable, x 

## Ask for degree and the matrix factorizations X and Y, then the hom 
## equiv classes, F and G 

n = input("Enter the degree of the singularity:") 

X1 = input("Enter the matrix factorization X, separated by a comma:") 

X = matrix([[0,X1[0]], [x1[1] ,0]]) 

Y1 = input("Enter the matrix factorization Y, separated by a comma:") 

Y = matrix([[0,¥Y1[0]], [Y1[1] ,0]]) 

Xn = matrix([[0,xX1[1]], [x1[0] ,0]]) 

Fi = input("Enter the equivalence class morphisms alpha and beta for F, 
separated by a comma:") 

F = matrix([[F1[0],0],(0,F1[1]]]) 

Gi = input("Enter the equivalence class morphisms alpha and beta for G, 
separated by a comma:") 


G = matrix({{0,G1[0]], [G1[1] ,0]]) 
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## Composes F and G into one matrix, called F, to then multiply 


F = F*G 


## Find the differential of Q associated to Y and multiply, then take 
## trace of new matrix 

dY = Y.derivative(x) 

mult = Fx*dY 


trace = mult(0,0] + mult[1,1] 


## Compute the derivative of singularity to find the denominator of 
## integral 

sing = x**n 

der = sing.derivative(x) 


divide = trace/der 


## Find the Taylor series expansion for the new polynomial and find 
## the coefficient of -1 degree term for final residue 

series = taylor(divide,x,0,-1) 

set = series.operands() 

1 = len(set) 

stset = str(set) k = len(stset) 


star = stset.find(’/x’) 


if star != -1: 
while star != -1: 
if 1 == 2: 


star = stset.find(’ ’) 


stset = stset[star+1:k-1] 
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break 
stset = stset[1:star] 
star = stset.find(’ ’) 
fin = SR(stset) 
print (fin) 
else: 


print") 
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